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I. INTRODUCTION 


A. BACKGROUND 


The art and science of marine positioning date back to antiquity, 
and have progressed from black magic, astronomical positioning, radar, 
Loran C and Omega to artificial satellite systems (Bresslau, 1980, pp. 
A35). The accuracy of land-based electronic systems varies: e.g., 
the Decca Pulse-8 system has an accuracy of + 50 m at 200 nmi; 
Loran-C has an accuracy of + 900 m at 2000 km (Laurila, 19/6, pp. 
348,429); Omega has an accuracy of +t 7.2 km globally (Creamer et al., 
1985,pp. 269); while the Navy Navigational Satellite System has a 30-m 
accuracy with precise ephemerides for a single pass. In general the 
accuracies of electronic positioning systems are degraded by land 
propagation distortion (Laurila, 1976, pp. 502). 

Astronomical navigation methods have serious observational limita- 
tions, i.e, the measurement depends on a plumb line and it is weather- 
sensitive, although a global datum is possible. The concept of 
positioning by an artificial satellite system began with the use of the 
Baker-Nunn camera for direction determination. This formed the foun- 
dation of real 3-dimensional positioning (Laurila, 1976, pp. 383). 

The current operational navigation satellites are part of the U.S 
Navy Navigation Satellite System, usually called TRANSIT. The system 
was conceived in 1958 at the Applied Physics Laboratory of the Johns 
Hopkins University and became operational in January 1964 and commer- 
cially available in 1967. The technical approach was to use the Doppler 
shift to provide intermittent position fix updates (Laurila, 1976, pp. 
463). The TRANSIT system consists of five operational satellites in 
polar orbits with altitudes of 1100 km, periods of 105 min, and frequen- 
cies of 150 MHz and 400 MHz. The satellites continuously broadcast 
predicted orbital coordinates along with timing signals, and are checked 
by four monitoring stations and two injection stations (Laurila, 1976, 
pp. 464). 


In addition, the Department of Defense has been developing the 
Navigation Satellite Timing and Ranging Global Positioning System 
(NAVSTAR GPS) under the Joint Program Office (JPO) since 1973 as a 
merger of Navy's Time and Navigation (TIMATION), with Air Force 
Project 621B (Bo0issenin, 1985, pp. 1). The basic method of 
GPS/TRANSIT position fixing differs from conventional celestial naviga- 
tion or visual satellite observation in that range, range-rate and angle 
are measured, since electronic systems do better in ranging than in 
angle measurement (Mueller, 1964, pp. 235-268; Easton, 1980, pp. 14). 
The TIMATION is a program to advance the development of high 
stability oscillators, time transfer, and two-dimensional navigation; 
Project 621b is to develop a highly accurate 3-dimensional navigation 
system (Boissenin, 1985, pp. 1). 

The development approach for GPS selected by the Defense System 
Acquisition Review Council (DSARC) was a step-wise development and 
test program, leading in successive phases to an operational GPS. Each 
phase was designed to build and expand on the previous phase in an 
integrated and cohesive manner. The decision of DSARC | in December 
1973 was to proceed with a phase | Concept and Validation to concen- 
trate on validation of design concepts through Development Test and 
Evaluation (DT&E) of user equipment. The decision of DSARC II in 
June 1979 approved proceeding into full scale engineering development 
to complete the DT&E and Initial Operational Test and Evaluation 
(lIOT&E) of user equipment. A DSARC III decision is planned for — 
October 1985 to obtain Approval for Full Production (AFP) of the user 
equipment (Boissenin, 1985, pp. 2). 


B. GPS POSITION AND NAVIGATION TECHNIQUE 


The Global Positioning System (GPS) is a space-based radio posi- 
tioning navigation system that is designed to provide continuous, highly 
accurate, three-dimensional position data to within a 16-m spherical 
error probability (SEP), velocity to within 0.05-0.10 m/s and system 


time to within 55 ns to suitably equipped users anywhere within 500 nmi 


of the earth. The position solution is based on World Geodetic System 
(WGS) 1972 ellipsoids in the Earth-Centered, Earth-Fixed (ECEF) coor- 
dinates. Accurate position and velocity determination will permit GPS to 
supplant less accurate systems such as Loran C, Omega, and Decca 
Pulse-8. GPS user equipment receives L-band signals from four satel- 
lites to measure pseudoranges and range-rates to solve for positions in 


three dimensions, velocity and time (Boissenin, 1985, pp. 2-3). 


C. GPS CONFIGURATION 


The GPS consists of space segments, control segments and user 


segments: 


1. §pgce Segment 


The space segment is composed of a constellation of satellites 
that continuously transmit signals to all users. When fully operational 
the GPS satellite constellation will deploy 18 satellites (plus 3 spares) in 
circular orbits at altitudes of 10900 nmi, inclined 45-65 deg, and having 
approximately 12-hr periods. The spacing of the satellites in their 
orbits will be arranged such that a minimum of four satellites will be 
available to a user anywhere on the globe (Boissenin, 1985, pp. 2; 
Jorgensen, 1984, pp. 1). Operational GPS satellites have projected 
mean mission durations of six years and design lives of seven and a 
half years (Senus and Heuerman, 1983, pp. 73). 

The unperturbed orbit of satellites may be given by six 
Keplerian elements, namely the semi-major axis (a), eccentricity (e), 
mean anomaly (M) or time of perigee passage (T), right ascension of 
ascending node (2), inclination of orbital plane with respect to the 
equatorial plane (!), and argument of perigee (w). The first three 
parameters determine the position of the satellite in its orbital plane, 
and the last three parameters determine the position of the orbital plane 
in the celestial coordinate system (Torge, 1980, pp. 115-119; Mueller, 
1964, pp. 147-157). 


tt 


The actual orbit of a satellite departs from the Keplerian 
model due to the effect of various disturbing gravitational forces and 
non-gravitational forces. The orbit can be viewed as oscillating ellipses 
which are given at each instant by the current orbital elements. The 
main gravitational forces are due to the earth, moon and sun. The 
earths gravitational force depends on satellite position and time due to 
irregularities in the earth's density. The non-gravitational forces are 
due to solar radiation pressure, magnetic effects, atmospheric drag, 
and relativity. The actual orbit can be modeled by the initial six 
Keplerian elements plus time (t), ellipsoidal harmonic coefficients (zonal 
harmonic J, tesseral harmonics J and K), and their rates of change 
(Torge, 1980, pp. 116-117; Heiskanen and Moritz, 1967, pp. 342-345). 
The computation of satellite position by ellipsoidal harmonic coefficients 
needs sophisticated software in the receiver (MacDoran et al., 1984, 
Bion. 205). 

The actual orbit can be modeled also by the initial six 
Keplerian elements, secular drift terms and harmonic coefficients to 
reduce computation. Such a model, however, requires continuous 
updating. As applied in the current GPS, corrections are uploaded to 
satellites every 24 hours by the master station, and these are used to 
make hourly changes on board the satellite (MacDoran et al., 1984, pp. 
59)5 

GPS satellites simultaneously transmit navigation information on 
two radio frequencies, L1 and L2. The L1 signal is modulated with C/A 
code (C/A = coarse-acquisition) and a P-code (P = precission), and the 
L2 signal is modulated with the P-code only. The functions of the 
codes are to identify the particular satellite and to measure phase 
shifts; both are necessary to measure pseudoranges. Both the L1 and 
L2 signals are also continuously modulated with the navigation data bit 
stream at 50 bps (Remondi, 1984, pp. 2; Milliken and Zoller, 1980, pp. 
6). Due to different gravitational potentials on the ground and at the 
satellite, and to velocities of users relative to the satellites, the satel- 
lites actually transmit a slightly lower frequency (Af= 4.45 x 107*° so 
that 10.22999999545 MHz is shifted to a nominal 10.23 MHz (Spilker, 
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1980, pp. 33). The navigation message contains information on the 
status of the satellite, the time synchronization information for transfer 
from C/A- to P-code, the parameter for computation of the clock 
correction, the ephemerides for computation of satellite position, the 
parameter for propagation delay correction, and the almanac for 
computing Geometric Dilution of Precession (GDOP)/satellite selection. 
The navigation message is intended to keep the satellite ephemerides 
and timing system as accurate as possible (Milliken and Zoller, 1980, 
pp. 7; Van Dierendonck et al., 1980, pp. 576). The contribution of 
satellite segment errors was monitored during the GPS development 
phase, and the recorded bias errors were found to be for X, -2 to -4 
fiptoGey, -o mm: and ftonmeZz, +3 m (JPO, 1984, pp. 6&). 


2. Control Seqment 


The control segment tracks all satellites in orbit. the 
purpose of satellite tracking is to observe the position and motion of a 
satellite as a function of time with a precision sufficient to refine its 
equation of motion to permit predicting future position with at least the 
accuracy achievable by direct observation. The operational control 
segment consists a master control station, a ground control station, 
monitor stations and an alternate control station. The monitor stations 
collect atmospheric and satellite tracking data from all satellites in view, 
which are encrypted and transmitted to the master control station, the 
ground control station and the alternate control station (Boissenin, 
1985, pp. 2). The master control station uses one week's data from all 
monitor stations to control the satellites’ health, to generate new satel- 
lite ephemerides and clock offset updates, and then transmit these 
updates to the constellation of satellites every 24 hours (Milliken and 
Zoller, 1980, pp. 69). 


3. User Segment 


The fundamental capabilities of the GPS user segment are 
satellite selection, signal acquisition, tracking and measurement, data 


recovery, and data processing for the real time user. The hardware 


if: 


component is comprised of an antenna, a receiver, and input/output 
devices (Glazer, 1980, pp. 81). 

If a GPS receiver has a very accurate clock synchronized 
precisely with the satellite clock, it can measure the transit time of the 
signal between the satellite and the receiver and multiply this by the 
speed of light to get the range from receiver to satellite. Due to 
expense and the impracticality of maintaining a synchronized clock at 
both the satellite and the receiver, only a crystal clock is used in the 
receiver, so that the measured signal transit time is subject to receiver 
clock error. This uncorrected distance is called the pseudorange. By 
measuring four pseudoranges to different satellites, the user's position 
can be calculated in the ECEF World Geodetic System 1972, which may 
then be transformed to another datum (Milliken and Zoller, 1980, pp. 
3-4). ; 

As the satellite signals vary, one can eceleit each signal by 
various means to get a position. At present there are two methods: 
The first uses the code signal to measure pseudorange and pseudorange 
rate; and the second uses the uncoded transmission signal as a source 
for very-long-base-line interferometry (VLBI). Boissenin (1985), 
Buennagel et al., (1984), Bock et al., (1984), Borel (1980), MacDoran 
(1984), JPO (1984), and Glazer (1980) provide detailed information on 
the uncoded GPS receivers. 

The Texas Instruments T1-4100 GPS receiver uses both 
C/A-code and P-code systems. It is a single-channel, two-frequency, 
digitized multiplexing receiver, capable of tracking simultaneously up to 
four satellites with continuous L1 and L2 frequency measurement. 
Lock-on for each signal is for 5 ms, for a total of 40 ms to complete 
one sequence (Johnson et al., 1984, pp. 61). The receiver clock 
permits connection to an external atomic clock so that the minimum 
number of satellites needed is reduced to three. The standard receiver 
unit consists of an antenna and pre-amplifier, a control display unit 
and a dual cassette recorder. The antenna assembly consists of a 
conical omni-directional antenna and the pre-amplifier. The Control 
Display Unit is a hand-held input/output device with alphanumeric 
displays (T1-4100, 1983, pp. 1-4 to 1-7). 
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Data acquired with a Tl-4100 receiver were used in this study 
which involved processing pseudorange data into position data and 
comparing GPS-determined positions with Trisponder-determined position 
data. The field data were acquired by the U.S Naval Oceanographic 
Office (NAVOCEANO) in Mississippi Sound. The development of the 
computer programs was very slow due to a large data set that was very 
difficult to handle; however, a mathematical model was developed. 
Calculations were actually made at the Naval Surface Weapon Center 
(NSWC) using NSWC computer programs written by S.L. Meyerhoff 
because the NSWC computer programs were found to be incompatible 
with the NPS IBM-3033; nor was it possible to get the NSWC computer 
program to run on the Fleet Numerical Oceanography Center's 
Cyber-1/75. The NSWC least squares adjustment program was more 
sophisticated than the one | developed in that it included a combination 
of bias Doppler data, bias frequency of the satellite's clock, and bias 


frequency of the receiver's clock. 
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ll. MATHEMATICAL MODEL USED IN THE POINT POSITIONING MODE 





The mathematical model may be divided into two parts, namely 
determination of coordinates in WGS 1972, and transformation to geodetic 
coordinates in some other datum. For this evaluation, the comparison 
between launch positions determined by Trisponders and by GPS is 
done in UTM, so a transformation of coordinates is performed for the 
shore station. 

To get good accuracy the relative geometry of the satellites and 
the ground receiver must be good. This geometry is calculated by the 
Geometric Dilution of Position (GDOP) method. The unknown position of 
the receiver's antenna (X, Y, Z coordinates) and receiver's clock offset 
are functions of four satellite positions and four corrected pseudo- 
ranges. The satellites’ positions are calculated from ephemerides. The 
receiver measures four observed pseudoranges to different satellites 
which must be corrected for atmospheric delay, relativity, earth rota- 
tion and satellite clock offset. The four unknown parameters of the 
receiver are calculated either by the conventional iterative least squares 
method or the Kalman filter method. Meyerhoff (1985), has developed 
computer programs to calculate either by least squares adjustment or 
Kalman filter methods. The formulas used for the computations in the 
computer program will be mentioned only here; their derivations can be 


found in the references. 


Ap SA)VEEEIVE SELEC MION 


If more than four satellites are above the user's horizon, the 
receiver can select the best geometry to get the best position, other- 
wise Satellite selection is made on the basis of the satellites health and 
the ages of the data messages. The geometric dilution of position 
(GDOP) is a measurement of the best geometry, which is a function of 


the user's and the satellite's position. The satellite's coordinates are 
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computed from the almanac. The parameters of the almanac are simpli- 
fied ephemerides parameters and are stored in the receiver's memory or 
can be entered by the operator; the formula used is similar to that 
developed in Subchapter C below. The computation of GDOP is 


performed in the receiver by the following formula: 


GDOP = [Trace (A A)-']*%? (1) 
COS dy COS ay Cos a7 | 
ie COS ay COS ay cos a% | 
cos a, COS ay Cos az 1 
coS a, Cos a, COsmma, i 
k a . — 
cos ay = direction cosines of the angles between the range to 


satellite number k and the x axis. 


An acceptable range is 2 < GDOP < 4 (Jorgensen, 1984, pp. 9-10). 


B. PSEUDORANGE 


The pseudorange corrections and satellites clock offset are calcu- 
lated and applied to the observed pseudoranges to get corrected pseu- 
doranges (Milliken & Zoller, 1980, pp. 4; Spilker, 1980, pp. 24). 


Ropbs= Rru* clt * (ot, - dt, )] 22) 
Rene = pseudo-range observed to SV-_ [ml] 
i ae = true range user to SV [m] SY ape’ 
ba = time correction consist of atmospheric, 


relativity, and earth rotation  [s] 


Site, = SV clock offset from GPS time [s] 
ota. = user clock offset from GPS time  [s] 
C = velocity of light [m/s] 
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The computation of satellite clock offset from GPS time is done by 
a polynomial model (Equation 2.4), given by Meyerhoff (1985) and Van 
Dierendonck et al. (1980). The model's parameters and the time refer- 
ence of the parameters are given in the navigation messages. The 
receiver's clock offset from the GPS time, is included in the four 


unknown receivers’ parameters. 


5tsy = clock offset correction of SV [s] 
teag = GPS time tag of measurement [fs] 
tee = the space vehicle code phase time at message 


transmission time fs] 


4 
= = ss = v4 
otey 997 21 *tag be) Y “Sas “Gey (2.4) 
a, = parameter of clock satellite bias  [s] 
a] = parameter of clock satellite drift [s/s] 
a, = parameter of clock satellite aging [s/s?] 
ae = reference time of clock drift at GPS system time  [s] 


The pseudorange correction consists of corrections due to atmos- 
pheric, relativity and earth rotation effects. Formulas for these are 
given in this section. The atmospheric correction may be separated 
into an ionospheric correction and a tropospheric correction (Equation 
2.5). The ionospheric correction depends on the electron contents of 
the atmosphere. The best ionospheric correction formula uses observa- 
tions of pseudorange at two frequencies Ll and L2. Equation 2.6 is 
used in the computation. The ionospheric correction may be calculated 
by an ionospheric model included in the satellite navigation messages; 
this polynomial model is less accurate than Equation 2.6, and was not 
used in this evaluation. If the P-code is not available, however, 


Equation 2.6 can't be used. 
fatm= tion® C trop (2.5) 
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The ionospheric correction for pseudoranges is given by Equation 
2.6 (Spilker, 1980, pp. 25), and the ionospheric correction for pseudo- 
range rates by Equation 2.7 (Meyerhoff, 1985, pp. 20). 


tea |. Josomor edLl (2.6) 

te = jonospheric correction  [s] 

StedL1 = difference between propagation at L1 and L2 [s] 
(a (8 ,/ Q- t/ Q,) (2.7) 

wlQhHY Qo)= - 1) 

tenia = ionospheric Doppler correction  [s] 

G = L1 Doppler count 

Q, = L1 frequency multiplier = 154 

g, = L2 Doppler count 

Q, = L2 frequency multiplier = 120 

yp = Nominal satellite frequency = 10,230,000 Hz 


The tropospheric correction is computed by the Chow model 
tropospheric correction given in Equation 2.8 (Meyerhoff, 1985, pp. 
22-24). 


— : 0.0276 Pp . (258) 
[sin 6 (tan 6 + 0.0045) + 0.00143] 
; A Omi Ese. 70s 10° “p 
[T? sin 6 (tan 6 +0.017) + 0.00035] 
E =e SOO el 10 
P = surface pressure [mbar] 
i = temperature [°K] 
0 = temperature lapse rate set at .006 [deg/m] 
H = percent relative humidity 
) = elevation angle of satellite [deg] 


The relativity corrections eae, due to the eccentricity of the 
satellite's orbit, are computed by Equation 2.9 (Meyerhoff, 1985, pp. 
27-28): 


Ile, 


t =e es sine wae (2.9) 


Cc 
Caan = relativity correction  [s] 
= eccentricity of satellite orbit 

E = eccentric anomaly (it could be solved 

iteratively) 
u = WGS-72 value of the earth's universal gravitational 

constant = 3.986008 10'* [m?/s?] 
A = semi major axis of orbit  [m] 
E = speed of light in vacuum [m] 


The earth rotation corrections (t, 54) due to the motion of the 
earth during the signal's propagation, are computed by Equation 2.10 
(Meyerhoff, 1985, pp. 25-26]: 


rote Ge C Ysat% ~ XsatYo) (2.10) 


Qe = WGS-72 value of the earth's rotational rate = 
0.00007292115147 [rad/s] 

Xsat Y%eat = ECEF of satellite coordinates, WGS-72 [m] 

yO) Es 10) <!5 coordinates of the station, WGS-72 [m] 


C. SATELLITE POSITION CALCULATION 


To compute a ship's position, the satellite position must be 
computed from the ephemerides transmitted in the navigation messages. 
As the satellite ephemerides parameters are Keplerian elements in 
secular drift terms and harmonic coefficients only, the receiver does not 
need sophisticated software to compute the satellite position. It must 
be noted, however, that the parameters change every hour. The 
computed mean motion was solved by Kepler's third law. The Keplerian 
parameters (mean anomaly, argument of latitude, radius of orbital 
plane, and longitude of ascending node) are affected by secular drift 
terms and harmonic coefficients, which must be corrected from their 


value at reference time by their drift. The drift of each parameter is a 


function of the time from epoch and secular drift plus harmonic correc- 
tions. A sample computation of satellite position may be found in 
Mueller (1964). 


1. Conversion of Receiver Time to GPS Time 


The GPS position is marked by the receiver clock; however, 
the reference time for satellite position computation is GPS system time. 
The receiver time tag must be converted to the GPS system time by the 


following step. 


Chae Veo SO 7 teran (2.11) 


t = GPS system time at time of transmission; GPS time of 


reception corrected for transit time  [s] 


fees = receiver time tag at receiver clock [s] 
dt, = receiver time offset from GPS system time originally set 


to zero, it is included in the unknown parameter [s] 
teran = transit time of signal from satellite to user [s] 


The time from epoch is related to the ephemerides reference time by: 


Ge tear tee (2.12) 


Fata 
1 


time from epoch, if greater than 302,400 subtract 
604,800; if less than -302,400 add 604,800 [s] 


tos = ephemerides reference time at GPS system time [s] 


2. Computation of Orbital Plane Parameters 





The orbital plane parameter is computed by Kepler's third law 
(Equation 2.13) and Kepler's first law (Equation 2.16). The WGS-72 


parameters for the computation are: 


m1 = WGS-72 value of the earth's universal gravitational 
constant = 3.986008 10'* [m?/s?] 
Q, = WGS-72 value of the earth's rotational rate = 


0.00007292115147. [rad/s] 
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The mean motion is affected primarily by secular drift terms; (Van 


Dierendonck, 1980, pp. 67) it is corrected as follows: 


5 = eee (2. 138 
A272 
No = computed mean motion [rad/s] 
A = semi major axis of orbit [m] 


A ae on (2.14) 
n = corrected mean motion [rad/s] 

dn = mean motion difference from computed value 

Mp NS ears (2.13% 
M = mean anomaly _ [rad] 
Mo = mean anomaly at reference time [rad] 


The eccentric anomaly was computed by: 


M =E +esineE (27 1G 


mM 
i 


eccentric anomaly (it could be solved 
iteratively ) 


~ = eccentricity of satellite orbit 
3. Computation of Orbital Plane Parameters in ECEF 


The orbital plane parameters in ECEF are affected by secular 
drift terms and harmonic coefficients; these effects were calculated as 


corrections to the parameters as follows: 


cos v = (cos E - e )/(1 - e cos E } (2.185 
sin v = [sin E (1 = e* )]7 Ol] e cosme) (2 1ah 
v = true anomaly [rad] 
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Cic 


r 


=yt wy (2.19) 


argument of latitude [rad] 


argument of perigee [rad] 
= C cos 2¢ + C,,sin 2¢ (2.20) 


argument of latitude correction [rad] 


amplitude of the cosine harmonic correction to the 


argument of latitude [rad] 


amplitude of the sine harmonic correction to the argu- 


ment of latitude [rad] 


ome eCOsuzor @csin 29 C221) 


radius correction [m] 


amplitude of the cosine harmonic correction term the 
orbit radius [m] 


amplitude of the sine harmonic correction term the orbit 


radius [m] 
Cnc Ccm Om aG , cin) 2¢ C272) 


= correction to inclination [rad] 

= amplitude of the cosine harmonic correction term the 
angle of inclination [rad] 

= amplitude of the sine harmonic correction term the 


angle of inclination [rad] 
= corrected argument of latitude [rad] 


=e ile] e cos E)*743* 6r (24) 


r = corrected [m] 

Pee i (2.25) 
| = corrected inclination [rad] 
L = inclination angle at reference time [rad] 

2 =2 + (2-2) t- % t (2.26) 


corrected right ascension of ascending node _ [rad] 


right ascension at reference time [rad] 


rate of change of right ascension [rad/s] 


Q 
a 
S 
4. Computation of Satellite Position in ECEF 


Equations 2.27 and 2.28 are used to compute the position of 
the satellite in its orbital plane (Mueller, 1964, pp. 147-170; Van 
Dierendonck et al, 1980, pp. 67): 


X= Tacos. U (2.2% 
Ye =r sin U (2.28) 
x = position in orbital plane  [m] 
NA = position in orbital plane  [m] 


Equations 2.29 and 2.30 are used to compute the position of the satellite 
ea a a al 


OT X* cos 2 - Y* cos | Sines (2.29) 
Memes <x sin 2% Y > “cosmicocme (2.30) 
Zee Na asiinel (2.318 


Xe atin ee 5 ECEF of satellite coordinate, WGS-72 [ml] 


sat sa 
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D. LEAST-SQUARES ADJUSTMENT (ITERATIVE METHOD) 


The corrected pseudorange between an unknown station P and the 


satellites is given by Jorgensen (1984): 


Roor = Reru * © (Sty - Stsv) (2.32) 
where: 
eae = corrected pseudorange from Satellite Vehicle (SV) 
number k to the unknown station P 
Bera = true pseudorange from Satellite Vehicle (SV) 
number k to the unknown station P 
D = VOX Bare. < neon 2° CZ 2 ee? (233) 
D = the preliminary distance between SV to the 


unknown station P. 
XsatYsat Zsat = the coordinates of SV [m] 
Xo, Yo, Zo = the preliminary coordinates of the unknown station 
P= [m] 
CG = the velocity of light in vacuum [m/s] 
The corrected pseudorange may be linerized by a Taylor series 


expansion as: 


Roor= Rea SRI AX*SRI AY*SR] AZ*ERI At 
5X] 5Y| 5Z| at] 


where the Jacobian of R calculated for 6X, 5Y, 6Z, and 6dt are 


B= 5R 5R 5R 5R 
5X 5Y §Z dt 
or 
Ng SEE G Moy £,~4sat c 
D 7) S) 


R = calculated pseudorange by preliminary coordinates of 
P from SV to the unknown station P- [m] 
The least squares adjustment may be done by the iterative method. 


The observation equation matrix: 


AB + v = f (2.34) 
f= Roar ~ Reor a 
where 
B = the Jacobian matrix, n by 4. 
A = the correction of the unknown coordinates matrix, 4 by 


1 matrix; Consists of 6%- AY | AZ =e anaes 


V = the residual observation matrix, n by 1] 
f = the difference between range corrected and range calcu- 
lated by a priori coordinates matrix, n by 1 
W = the weight matrix, n by n 
w = 1/0; o is the variance of the observations 
n = number of observations, minimum 4 observations 
7 ae calculated pseudorange matrix by preliminary coordinates 
of P from SV to the unknown station P 
Reor = corrected pseudorange matrix SV from the unKnown 


station P 
Minimizing Equation 2.34 as is done by Mikhail and Ackerman 
(1974) 


N = B'wB (2,365 
- Ak 
t = B wf (2-3 
A= N7 t= (B-WB)7 (B we) (2.38) 


Xp, \p, Zp = the corrected coordinates of the unknown station P  [m] 
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Xp = % * AM) (39) 
mee) Yu) | AC) (2.40) 
Zee * (ACS) (2.41) 
Ot yaaa) (2.42) 
Xo, Yo, Zo = the preliminary coordinates of the unknown station P 

[m] 

A = estimator of reference variance matrix 

ACi) = element of matrix A [m] 


E. KALMAN FILTER 


Meyerhoff (1985) gave a method for the computation of dynamic 
positions by means of an eight-state Kalman filter after pseudoranges 
have been corrected and smoothing has been applied to them by Doppler 
count. The computation of position consists of computation of corrected 
pseudoranges, satellite's positions, variance of pseudoranges, Jacobian 
matrix B, Kalman gain (K), updated position, and the updated covari- 
ance matrix (P). [For details of the Kalman filter see Meyerhoff 
(1985), Bierman (1982), and Brown (1983)]. The variances of the 
pseudoranges are computed from the variance pseudoranges of L1 and 


L2. The variances of the pseudoranges of L1 and L2 are 


val = Ko Bs K; [ Br. n ] [ Che. * Bp] (2.43) 
10°L1 10°L1 
ll = variance of pseudorange data for L1 or L2, respectively 


at timetag  [s?*] 


Ki&K2 = variance bias factors’ [s?] 

BF = code loop noise bandwidth 

n = number of trackers 

Gi = 0.1 of carrier power to noise density of L1 
C1 = 0.1 of carrier power to noise density of L2 
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By = code predetection bandwidth 


The variance of pseudorange is computed by: 


Vosd = [Vui* M1 * M2 7 £0Q,/Q, 3? -1)? (2.44) 
Vosd = cariance of pseudorange '  {f[s?] 
E(k, kK) = Nene for tracker k, respectively 
K = P(-) B'(B P(-) BE + 5)7? (2.45) 
K = Kalman gain matrix 
r = variance of pseudorange matrix 
B = Jacobian matrix 
P(-) = The covariance matrix before updated 
of this data point, the initial value is the initial vari- 
ance of X , Y , Z, 6&t ; of diagonal are zero 
PO)S = Pl) hae (2.46) 
P(+) = The updated covariance matrix 
Rome. > Robs ®eai “© * cor (2.47) 
Reon = pseudorange range residuals [ml] 
Robs = pseudorange observed - [m] 
eT = pseudorange calculated [m] 
eee = correction of pseudorange’- [m] 
c = speed of light [m/s] 
rE ree CD) ke (2.48) 
Y (eae) ae (2.49) 
X(-) = X coordinate before update [ml] 
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Gas Ae IK) Rome 


5t(+) = 


vi) 
Zi(=) 
6t(-) 
K Ci) 
no G) 
rece) 
Za) 
1, 


mej) 1) a ie 


Y coordinate before update  [m] 

Z coordinate before update [ml] 
receiver clock offset before update  [s] 
row i of matrix K 

The updated X coordinate [ml] 

The updated Y coordinate [ml] 

The updated Z coordinate [m] 


The updated receiver clock offset [s] 


Ia, 


(2.50) 


(eo 1) 


ltl. THE TIl-4100 TEST IN MISSISSIPPI SOUND 





A. FIELD OPERATION 


The U.S Naval Oceanographic Office (NAVOCEANO) conducted a 
test on Mississippi Sound in May 1984 to evaluate the performance of 
the TIl-4100 on a marine platform (small launch) and to determine 
whether the T1I-4100 is adequate for hydrographic surveying. The 
launch position was determined by a GPS receiver and four Del Norte 
Trisponders on shore stations; two of the Trisponders were located 
above second-order conventional North American Datum-192/7 (NAD-27) 
bench marks, while the other two Trisponder positions were determined 
by means of the TRANSIT SYSTEM (Dupont and Dunn, 1984, pp. 11-49 
and 11-50). One TI-4100 was on the launch and one was at a third 
NAD-27 bench mark. Three first-order bench marks were observed by 
a TRANSIT receiver to determine the local datum shift from NAD-27 to 
WGS-72 (Dunn, 1985). 

Data for 1984 Julian days 136 and 138 collected during the 
NAVOCEANO test were arranged in standard GPS data exchange format 
(Scott and Peters, 1983, pp. 3-23) and sent to NPS for study. The 
GPS shore station data and Trisponder data were stored on a single 
data tape. The GPS shore station data were acquired every 12 
seconds, and the Trisponder data were acquired every second. The 
launch receiver's raw GPS data (consisting of pseudoranges, 
pseudorange-rates, satellites ephemerides; timetag and atmospheric 
data) were acquired and recorded every second; the launch's real-time 
positions as determined by GPS (X, Y, Z and 6t) were acquired every 
6 seconds; and both the raw data and real-time positions were stored 


on three additional tapes. 
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B. DATA PROCESSING 


NSWC has four computer programs for post-processing raw GPS 
data into the receiver's position and clock offset. The first program is 
called CONSTR; the second SOLTOM; the third KALMN2; and the fourth 
program DOALL. The first converts data from the standard data 
exchange format to a new format compatible with the second and third 
programs. The SOLTOM and KALMN2 programs convert the output data 
of CONSTR into ECEF coordinates by least square adjustment or by 
eight-state Kalman filtering, respectively. DOALL compares dynamic 
post-processed data with the Trisponder position data in UTM coordi- 
nates. The programs were written by Meyerhoff (1985) for the 
Computer Data Corporation Cyber-1/75 and in their entirety are quite 
large. | was not able to convert these programs for the NPS IBM-3033 
due to their complexity and size. | 

The data from GPS time 185000 s to 188600 s (At = 3600 s) on 1984 
Julian day 136, was converted from standard 9-track tape to a new 
format, so that it would be compatible for further processing by the 
least squares adjustment or Kalman filter computer program. | 

For the 3600-s interval for which shore station data were analyzed 
using the least squares adjustment program, the positions obtained by 
post-processing did not significantly differ from real time positions. 

For the 3600-s interval for which GPS and Trisponder launch posi- 
tions were analyzed, a total of 31/1-s data points from both positioning 
systems were available. The launch positions determined by GPS 
consisted of 2/40-s of data points based on four satellites and 431-s of 
data points based on three satellites. The positions of the launch 
determined by GPS were post-processed using the Kalman filter program 
and compared with the launch positions determined by Trisponders. 
Since the time tags of the shore station's GPS positions differ from the 
launch's, further differential mode post-processing of this data was not 
done. 

The post-processing for GPS shore station positions and launch 
positions used broadcast ephemerides. Post-orbit ephemerides were not 


available to possibly improve the accuracy. 
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C. ACCURACY OF LAUNCH POSITION DETERMINED WITH 
TRISPONDERS 


The "true" position of the launch was computed by least squares 
adjustment from the two second-order shore stations (NAD-27) and the 
two TRANSIT stations in WGS-72 using the four ranges from the 
Trisponders. The two second-order shore stations (NAD-27) were 
transformed to WGS-/72 and converted to UTM. The two TRANSIT 
stations in WGS-72 were converted from ECEF to UTM. The total accu- 
racy of the launchs position determined by the Trisponders is + 7 m, a 
value which depends on the accuracy of the Trisponder locations and 
the accuracy of the Trisponders themselves; Table | shows the compu- 


tation. 


TABLE | 


THE TOTAL. ACCURACY OF THE LAUNCH'S POSITION 
DETERMINED BY TRISPONDERS 

The accuracy of transformation = 2°6 Mm 
The accuracy of Trisponder = 2 3m 
The accuracy of convertion to UTM = + 1m 
The accuracy of Transit =. = “Nem 
The rms accuracy of the 
launch position determined 
by Trisponder = + 6.86 m 


The accuracy of the launch's position determined by the 
Trisponders depends on the accuracy of the position of the two second- 
order Class 2, shore stations used as Trisponder sites locations. For a 
conventional geodetic station Bossler (1984) gives second-order accuracy 
as 1:20,000. The length of a side is generally not less than 5000 m, 
and the standard deviation of angle observation does not exceed 0.8". 
Thus, the position accuracy of the second-order shore station sites is 
approximately + 0.25 m. The accuracy of the transformation from 
NAD-27 to WGS-72 was determined by datum shift. The datum shifts 
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were determined locally by TRANSIT observations on three first order 
NAD-27 stations. The root-mean-square accuracy of the transformation 
may be estimated as approximately 5 to 10 m (Kumar, 1985). The accu- 
racy of conversion from WGS-/2 to UTM coordinates may be assumed + 
1.0 m (Kumar, 1985; U.S Army, 1973, pp. 8-9) 

The accuracy of the launch's position determined by the 
Trisponders depends also on the position accuracy of the other two 
shore stations determined by TRANSIT satellite observations. The 
accuracy of the TRANSIT stations may be assumed to be + 1.0 m, since 
observations from more than 25 passes were used, and precise ephemer- 
ides were applied (Bossler, 1984, pp. 3-10 and 3-11; Kumar, 1985). 

Because the Trisponder is a microwave positioning system, it is 
subject to multipath errors and range holes. In general the accuracy 
of a Trisponder is + 3 m (Umbach, 1981, pp. AD-2), although it might 
be larger in some environments. Because the Trisponder observations 
were made at night during satellite passes, the accuracy of the test 


data may be better than + 3 m. 


D. ACCURACY OF THE LAUNCH'S POSITION DETERMINED BY THE 
GPS RECEIVER 


The GPS post-processing position values were compared with the 
launch's position determined by the Trisponders, and the mean of the 
rms values was found to be less than + 11 m if four satellites were 
available (Table I!). Positions for the first three seconds were poor, 
with the maximum rms difference being 215 m. The mean of the rms 
values was found to be less than + 21 m if only three satellites were 
available (Table III). 

The signal quality of satellite number 8 was variable; this was 
indicated by bad pseudorange residuals, which caused an increasing of 
the northing, easting and rms differences. The signal quality of satel- 
lite number 8 tended to be cyclic, decreasing until it was lost, leaving 
only 3 satellites, and then increasing again, making four satellites 


available again. The signal quality is measured by the pseudorange 
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TABLE II 
DIFFERENCE BETWEEN a AND IR 


im SPONDER-DETERMINED 
POSITIONS WIT IN V 


| 
LEMES 


FOUR SAT IEW 
Mean Maximum Minimum Sigma 
rms 10267 7 250035 0.422 125 735 
Easting -/.058 12e7on =202 622 10.653 
Northing 2.136 25.445 -19.429 6.854 


residuals, which are considered bad when exceeding 20 m. A total of 
865 available data points indicated bad pseudorange residuals, which 
were the result of bad Doppler counts. The number of bad Doppler 
counts were: 1 from satellite no. 6, 1 from satellite no. 3, 97 from 
satellite no. 11, and /66 from satellite no. 8. 

The presence of the four "good" satellites and four good pseudo- 


ranges will significantly improve the accuracy. 


TABLE 111 
DIFFERENCE BETWEEN GP AND TRISPONDER-DETERMINED 
POSITIONS WITH THREE Sa TELLITES IN VIEW 
Mean Maximum Minimum Sigma 
rms 20). 906 133762) 0.765 Bigeo3d 
Easting =19. 920 Te lez -125.801 2051938 
Northing -1.614 324256; -45.040 ze 


The accuracy of the launch's position determined with the 
- Trisponders is better than the launch's position determined with GPS. 
Although 29 percent of the data have similar accuracy, further studies 
to separate good data may be performed. The T1-4100 does not, 
however, have the distance limitation of the Trisponders: It is limited 
only by the present coverage of GPS satellites. 

The real-time dynamic positions were compared with position deter- 
mined by post-processing. The mean position differences were AX= 80 


m, AY= 7m, and AZ= 13 m. 
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IV. CLUSIONS 


1. The overall accuracy of Trisponder-determined positions for this 
study was + / m. As the observed range increases, the accuracy of 
positions obtained from the Trisponders decreases. Beyond line-of- 
sight Trisponders cannot be used. 

2. The accuracy of ship positions determined at sea by a T1-4100 
GPS receiver is slightly lower than the position accuracy obtained with 
shore-based Trisponders, although 29 percent of the data analyzed have 
similar accuracy. The GPS receiver is more versatile than the 
Trisponder system in that GPS is not limited by line-of-sight. Post 
processing accuracy may be improved if post-orbit ephemerides data 
were available for the analysis. 

3. To improve the position accuracy by GPS receivers, further 
studies in a differential mode should be performed in the marine envi- 
ronment. | 

4. The main difficulty of this study was the complexity of the 
programs for the CDC computer. The conversion of the NSWC computer 
program to the NPS IBM-3033 is suggested for the further application at 
NPS. A collection of hydrographic subroutines in a common library at 


NPS would be very helpful in future hydrographic studies at NPS. 


SO 


AFP 
C/A-code 
Coe 
DSARC 


DT&E 
ECGEE 
FNOG 
GDOP 
GPS 
IBM 
IOTSE 
JPO 
MHz 
NAD-27 
NAVOCEANO 
NPS 
NSWC 
P-code 
rms 

SEP 

SV 
Timation 
USN 
UTM 
VLBI 
WGS -72 


APPENDIX A 


ABBREVIATIONS AND ACRONYMS 


Approval for Full Production 
Coarse and Acquisition code 
Computer Data Corporation 
Defense System Acquisition Review 
Council 

Development Test & Evaluation 
Earth Centered Earth Fixed 

Fleet Numerical Oceanography Center 
GEOMETRIC Dilution of Precision 
Global Positioning System 
International Business Machine 
Initial Operational Test and Evaluation 
Joint Program Office 

Mega Hertz 

North American Datum 1927 

U.S Naval Oceanographic Office 
Naval Postgraduate School 

Naval Support Weapon Center 
Precision code 

root mean square 

Spherical Error Probability 
Satellite Vehicle 

Time and Navigation 

United States Navy 

Universal Traverse Mercator 
Very-Long-Baseline-Interferometry 
World-Geodetic-System 19/72 
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